Finite-Size Scaling Study of the Surface and Bulk Critical Behavior 
in the Random-Bond 8-state Potts Model 
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The self-dual random-bond eight-state Potts model is studied numerically through large-scale Monte 
Carlo simulations using the Swendsen-Wang cluster flipping algorithm. We compute bulk and surface 
order parameters and susceptibilities and deduce the corresponding critical exponents at the random 
fixed point using standard finite-size scaling techniques. The scaling laws are suitably satisfied. We 
find that a belonging of the model to the 2D Ising model universality class can be conclusively ruled 
out, and the dimensions of the relevant bulk and surface scaling fields are found to take the values 
yh = 1.849, yt = 0.977, yi ls = 0.54, to be compared to their Ising values: 15/8, 1, and 1/2. 
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The understanding of the role played by impurities on 
the nature of phase transitions is of great importance, 
both from experimental and theoretical perspectives. It 
is a quite active field of research where the resort to 
large-scale Monte Carlo simulations is often necessary jl| . 
The effect of quenched bond randomness in a system 
which undergoes, in the homogeneous case, a second- 
order phase transition has been extensively studied. It is 
well understood since Harris proposed a relevance crite- 
rion for the fluctuating interactions || . Disorder appears 
to be a relevant perturbation when the specific heat ex- 
ponent a of the pure system is positive. Since in the 
two-dimensional Ising model (IM) a vanishes due to the 
logarithmic Onsager singularity, this model was carefully 
studied in the '80s [§. 

The analogous situation when the pure system ex- 
hibits a first-order transition was less well studied, in 
spite of the early work of Imry and Wortis who ar- 
gued that quenched disorder could induce a second- 
order phase transition [Q. This argument was then 
rigourously proved by Aizenman and Wehr, and Hui 
and Berker ||. In two dimensions, even an infinitesi- 
mal amount of quenched impurities changes the transi- 
tion into a second-order one. The first intensive Monte 
Carlo study of the effect of disorder at a first-order phase 
transition is due to Chen, Ferrenberg and Landau (CFL). 
These authors studied the q = 8-state two-dimensional 
Potts model, which, in the pure case, is known to ex- 
hibit a first-order phase transition when q > 4 ||. They 
definitively showed that the transition becomes second- 
order in the presence of bond randomness JtJ, and ob- 
tained the critical exponents from a finite-size scaling 
(FSS) study at the critical point of a self-dual disordered 
system [p[. Their results, together with other related 
works |lO] , |ll|] , suggested that any two-dimensional ran- 
dom system should belong to the 2D pure IM universality 
class |Q. In a recent paper, Cardy and Jacobsen (CJ) 
used a different approach |13| ], based on a transfer ma- 
trix formalism vM, to study random-bond Potts models 



for different values of q. Their estimation of the critical 
exponents leads to a continuous variation of /3b /v with 
q. This result is in accordance with previous theoretical 
calculations when q < 4 Jig] , and, in the randomness- 
induced second-order phase transition regime q > 4, it 
is quite different from the Ising value of ^, and in sharp 
disagreement with the Monte Carlo results of Ref. Q for 
g = 8. 

The surface properties of dilute magnetic systems paid 
less attention. Quite generally, the scaling laws involving 
surface and/or bulk exponents can be deduced from a ho- 
mogeneity assumption for both surface and bulk singular 
free energies e.g. 

f surf (t, h, h.) = b-^f sarf {b y % b^h, b v »° h s ). (1) 

All the standard critical exponents can be expressed in 
terms of the anomalous dimensions yi associated to the 
relevant scaling fields This makes their determina- 
tion of great importance in the case of random systems. 
The (1,1) surface of the Ising model on a square lattice 
has only recently been investigated through Monte Carlo 
simulations by Selke et al jD]]. The critical exponent j3\ 
of the surface magnetization was found to be very close 
to its value in the pure 2D IM. 

In this letter, we report a FSS study of the bulk and 
surface critical properties of the 8-state random-bond 
Potts model. Although this model has already been stud- 
ied by Monte Carlo simulations, our approach is the first 
investigation of the surface properties for a random sys- 
tem other than Ising-like. It leads furthermore to dif- 
ferent results in which concerns the bulk properties, and 
our aim is to bring some clear evidence to solve the dis- 
crepancy between the recent results of C J (Ref. ) and 
those of CFL (Ref. §). 

In the following, we consider the q = 8-state random- 
bond Potts model on the square lattice. The Hamilto- 
nian of the system with quenched random interactions is 
written 



1 



2 



CHRISTOPHE CHATELAIN AND BERTRAND BERCHE 



(2) 



where the spins take the values a — 1,2, ... ,q and the 
sum goes over nearest-neighbor pairs The cou- 

pling strengths are allowed to take two different values 
K\ = K and K2 = Kr with probabilities p and 1 — p 
respectively. The ratio K2/K1 is kept to the constant 
value r = 10 (a strong enough value in order to ensure 
that the critical behavior is no longer governed by the 
pure system fixed point) for all the simulations. If both 
couplings occur with the same probability, p — 0.5, the 
system is, on average, self-dual, and the critical point is 
exactly given by the critical line of the usual anisotropic 
model |l|Jl|: (e K ° - l){e K " r - 1) = q. At p = 0.5, 
we performed large-scale simulations of L x L lattices 
(10 < L < 96) with periodic boundary conditions in one 
direction (vertical direction) and free boundaries (FBC 
system) in the other direction. An equal number of cou- 
plings of both type is first distributed over all the bonds 
of the lattice. The couplings are then mixed randomly. 
This procedure ensures an exact probability p = 0.5, and 
avoids the fluctuations around this value which would 
result from the use of a random number generator to 
build the distributions of couplings. The multi-spin cod- 
ing and the Swendsen-Wang cluster flipping method |l^] 
were used, and the histogram technique allowed us to 
determine the behavior of the different quantities over a 
range of K ^0|. For each distribution of the couplings, 
between 2 x 10 5 (smaller lattice sizes) to 4 x 10 5 (larger 
lattice sizes) Monte Carlo steps per spin were performed. 
Although it is smaller than the calculations of Ref. 
this is always larger than 10 4 times the correlation time, 
and turns out to be sufficient in order to produce reliable 
thermal averages. On the other hand, around 30 disor- 
der realizations were performed in Ref. but, since the 
averages over randomness are still strongly fluctuating 
(Fig. |l|), we used 500 (10 < L < 32), 330 (40 < L < 64), 
and 250 (72 < L < 96) disorder realizations. (For smaller 
lattice sizes, the configurational fluctuations of data due 
to randomness are more pronounced, and a larger num- 
ber of configurations is needed). These values guarantee 
the same order of magnitude for the contributions to the 
statistical errors resulting from the thermal average and 
from the replica average. 

The translational invariance is restored in the vertical 
direction by averaging over the disorder realizations. The 
local order parameter for a given replica, written fij, is 
defined by the majority orientation of the spins at column 



g-i 



(3) 



Here, p ma x(j) = max a (p<?(j)), where p a {j) is the den- 
sity of spins in the state a at column j. The thermal 
average over the Monte Carlo iterations, written with 



brackets (...), is performed and the physical quanti- 
ties are then averaged over disorder configurations, for 
example m,j = [{fij}], where [. . .] denotes the replica 
average. The local surface susceptibility is given by 
X11 = KL \{pi) — (mi) 2 ] an d similar quantities are de- 
fined for the bulk. 
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FIG. 1. Fluctuations of the averages over the number of re- 
alizations of disorder. Average total magnetization for L = 56 
over up to 330 replicas. 

The first part of our analysis was to study the lo- 
cal magnetic surface properties. The local surface mag- 
netization mi = [{p\)\ is expected to follow the usual 
finite-size scaling behavior at the infinite lattice critical 
point: mi(K c ,L) ~ L~ Xl , where the critical dimension 
x i = Pi/v is deduced from the size dependence of mi, 
[(Mi 2 )], and [< Ml 4 >] (Fig.|). 
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FIG. 2. Log-log plot of [(/^i)] (circles), [(/i] 2 )] (squares), 

and [(mi 4 )] 1,/4 (diamonds) vs L. The slopes over the whole 
range of values of L are respectively —0.461, —0.456, and 
—0.443. The insert shows the effective surface exponent 
xi(L m in) defined in the text vs L^ in (deduced from mi) and 
its extrapolated value. 

A power-law fit of the curve between a given smaller 
size L m in and the maximal value L max = 96 defines an ef- 
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fective exponent x\(L m i n ). The smaller size is then can- 
celled from the data and the whole procedure is repeated 
until the three larger sizes only remain. The effective ex- 
ponent is plotted against £ m *„ (Insert in Fig. ||), and the 
critical exponent follows from the extrapolation at infi- 
nite size in the linear regime. Here the final estimation 
gives: 



Pi 



0.467 ±0.006, 



(4) 



where the uncertainty is the standard deviation. 

The behavior of the local surface susceptibility \u is 
more ambiguous, since xn seems to exhibit a power-law 
behavior with a very small exponent, but also fits with a 
logarithmic divergence as it is the case for the pure IM. 
From the behaviors of xn{K c ) and xTi ax (deduced from 
histogram reweighting) , we obtain Ju/p = 0.099 ±0.009 
(Fig. ||) . It corresponds to a fit which gives greater place 
to large sizes. On the other hand a logarithmic behavior 
seems also convincing (Fig. ||). Although this first anal- 
ysis does not allow any definitive conclusion, the scaling 
relation 2f3±/i> + ~fn/v — d — 1 is best satisfied, within 
error bars, by the power-law case and is furthermore sat- 
isfying with the value of x\ which rules out the 2D IM 
universality class. 
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FIG. 3. Local susceptibility ^ii(^c) at the critical point 
vs L and power-law fit (solid line) or logarithmic fit (dashed 
line). The insert shows the effective exponent ju/v. 

Once the local surface properties have been studied 
and since they strongly suggest that the model does not 
belong to the 2D IM universality class, a carefull analysis 
of bulk properties is needed to bring definitive conclu- 
sions. Bulk properties are furthermore needed in order 
to compute surface excess magnetization. For this pur- 
pose, we made simulations on L x L lattices with periodic 
boundary conditions in both directions (PBC system). 
The average quantities over the whole system lead to the 
critical exponents associated to the bulk magnetization 
and bulk susceptibility (Fig. ||). The determination of 



the slopes in the log-log plots, and of the correspond- 
ing standard deviations (of the order of 10%), indicates 
large fluctuations. We then turned back to the "effective 
exponent" technique presented above to have accurate 
estimations. One thus obtains 



— = 0.153 ±0.003, — 
v v 



1.701 ±0.008. 



(5) 



The first value is very different from the IM value and 
from the result of Ref. Q (0.126) and closer to the result 
of Ref. (l| (0.142). The second value satisfies to bet- 
ter than 0.7% the scaling law resulting from Rushbrooke 
and hyperscaling relations: 7t/f = d — 2f3b/v. The cor- 
relation length exponent v is deduced from the devia- 
tion of the effective critical coupling (at the maximum 
of the bulk susceptibility %™ aa: ) from its exact value, 
\K C (L) - K c \ ~ L~ x l v (Fig. f§. It leads to the corre- 
lation length exponent v = 1.023 ± 0.020 (Fig. |) [§2|. 
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FIG. 4. Bulk quantities: Maximum of the susceptibility 
(Xb lax , squares) and its value at K c (xb, circles); average mag- 
netization (mt, diamonds); deviation of the effective critical 
coupling from its exact value (triangles), and surface excess 
magnetization (m 3 , crosses). The corresponding exponents 
7b/j/, —f3b/i>, —1/v, and —f3 a /v are given in the text. 

Excess surface magnetization can be calculated by a 
comparison between the FBC and PBC systems: 



1 L 



1 



Tfih — rrij) ~ —L{rnb — m) 



(G) 



where fh is the average magnetization for the FBC sys- 
tem and mb for the PBC one. In Eq. |^, the approxima- 
tion symbol renders the possible difference between the 
majority spin orientation on a layer j and its value for 
the whole system. It produces a small difference between 
mean magnetization and mean profile for the FBC sys- 
tem but on average the two quantities should scale the 
same way. The corresponding exponent obtained from 
FSS takes the value: 
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— = -0.852 ±0.007. 
v 



(7) 
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It is in accordance with the expected value resulting from 
the scaling law (5 s /v = 0b/v — 1 = —0.847, albeit m s , 
given by a difference, is subject to strong fluctuations. 

In this letter, we reported the results of large-scale 
Monte Carlo simulations [E3| of the surface and bulk 
critical behaviors of a randomness-induced second-order 
phase transition in the 8-state Potts model. Concern- 
ing the bulk critical exponents, there is a clear discrep- 
ancy between our results and those of Ref. || which were 
very close to the pure IM values. The main difference 
between our procedure and these previous simulations 
is due to our use of a larger space of allowed coupling 
strengths. The self-dual condition p = 1/2 was indeed 
imposed by these authors in both directions. A possible 
weakening of randomness could result of this choice. We 
furthermore generated a number of disorder realisations 
10 times larger, which makes our results reliable. On the 
other hand, our value for (3b /v is slightly above Cardy and 
Jacobsen's result, while we used the same order parame- 
ter than CFL. The possible explanation suggested in CJ 
(Ref. |l3|]) for the disagreement with CFL (non-standard 
order parameter) has thus to be dismissed. 

We have also reported here the first extensive numeri- 
cal study of surface critical behavior in a randomness- 
induced second-order phase transition. While excess 
magnetization offers an alternate determination of the 
scaling dimension of the bulk magnetic field, local sur- 
face properties lead to the scaling dimension of a surface 
field which is also relevant. 

We can summerize all the results in a table of the 
anomalous dimensions of the relevant scaling fields (Ta- 
ble ||). The independent determinations of these values, 
very close together, give reliability to the results. The fi- 
nal estimations are the following: y t = 0.977, yt = 1.849, 
Vh. = 0.54. 

TABLE I. Scaling dimensions of the bulk and surface fields 
and of the temperature deduced from the values of the critical 
exponents. 
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Finally, one has to mention that we also computed 
profiles and correlations (details will be published else- 
where). The values of ry = 0.29 (correlations PBC sys- 
tem), and of the critical exponents difference x\ — Xf, — 
0.27 (profile close to the free surfaces, FBC system) lead 
to results which are slightly too small compared to the 
previous values of X\ and Xb- Surprinsingly, the bulk ex- 
ponent Xb — 0.145 is found to be very close to Cardy and 
Jacobsen's result which was deduced from the behavior 
of correlations as well, but within a strip geometry. 
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